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Abstract: We investigate the localization problem of matter fields within the 5D stand- 
ing wave braneworld. In this model the brane emits anisotropic waves into the bulk with 
different amplitudes along different spatial dimensions. We show that in the case of increas- 
ing warp factor there exist the pure gravitational localization of all kinds of quantum and 
classical particles on the brane. For classical particles the anisotropy of the background 
metric is hidden, brane fields exhibit standard Lorentz symmetry in spite of anisotropic 
nature of the primordial 5D metric. 

Keywords: Field Theories in Higher Dimensions, Large Extra Dimensions, p-branes 
ArXiv ePrint: 1204.2448 [hep-th] 



Contents 



1 Introduction 



1 



2 Background solution 



2 



3 Localization of bosons 



4 



4 Localization of fermions 



7 



5 Classical particles 



10 



6 The summary and discussions 



12 



1 Introduction 

Braneworld models involving large extra dimensions [1-4] have been very useful in ad- 
dressing several open questions in modern physics (for reviews see [5-8]). Most of these 
models were realized as time independent field configurations. However, there have ap- 
peared also several non-stationary braneworlds [9-12]. Here we consider the scenario with 
time-dependent metric coefficients recently proposed in [13, 14]. 

A key requirement for realizing the braneworld idea is that the matter fields be localized 
on the brane. For several reasons one would like to have a universal gravitational trapping 
mechanism for all fields. However, there are difficulties to realize such a mechanism with 
exponentially warped space-times. In the existing 5D models the spin and spin 2 fields 
can be localized on the brane with the decreasing warp factor [3, 4], while the spin 1/2 
fields can be localized with the increasing factor [15], and spin 1 fields are not localized 
at all [16]. In the case of 6D models it was found that spin 0, spin 1 and spin 2 fields 
are localized on the brane with the decreasing warp factor and spin 1/2 fields again are 
localized with the increasing factor [17]. There exist also 6D models with non-exponential 
warp factors providing the gravitational localization of all kinds of bulk fields on the brane 
[18-20], however, these models require an introduction of unnatural sources. 

In this paper we show the universal gravitational trapping of all kinds of matter fields 
within the 5D standing waves braneworld [13, 14]. We start (Sect. 2) from recalling the 
basic ingredients of the model. In Sect. 3 and 4 we demonstrate the existence of zero 
modes of all kinds of quantum particles on the brane. The Sect. 5 is devoted to the 
geodesic motions of classical particles. A short summary and discussion are presented in 
Sect. 6. 
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2 Background solution 



In this section we briefly review the 5D standing wave braneworld model [13, 14]. The 
braneworld is generated by gravity coupled to a massless ghost scalar field, which depends 
on time and propagates in the bulk. In order to avoid the well-known problems of stability 
which occur with ghost fields, the bulk scalar field does not couple to ordinary matter in 
our model [13, 14]. For other models with phantom fields see [21-25]. 
The action of the model has the form: 



S = / d'x^ 



^ {R-2A)+g^''dMcl>dN^ 



IGttG 



(2.1) 



where the capital Latin indexes refer to 5D space-time and G and A are the 5D Newton 
and cosmological constants, respectively. 
We use the metric ansatz: 

ds^ = e^^l^'l {dt^ - e'^dx^ - e'^dy^ - e'^^dz^) - dr^ , (2.2) 
where a is a curvature scalar, with the determinant: 

Vff = e4'^H. (2.3) 

The system of coupled Einstein and the scalar field equations for the ansatz (2.2) is 
self-consistent if: 



n(t,r) = y-7rG(/)(t,r) , A = 6a\ (2.4) 
and it has the standing wave solution [13, 14]: 

u{t,r) = sm{ujt)Z{r) , (2.5) 

where 

Z(r) = Ce-2a|r| J^^g-a|r|^ _ (2.6) 

Here C and oj are real constants and J2 is the second-order Bessel function of the first kind. 
Thus the model (2.2) describes the brane located at r = which possesses the anisotropic 
oscillations and sends a wave into the bulk (as in [26, 27]), i.e. the brane is warped along 
the spatial coordinates through the factors ~ e"*-*'*"-*, which depend on time t and the extra 
coordinate r. 

As pointed out in [13, 14], the ghost-like field (j){t, along with the metric oscillations 
u{t,r), must be unobservable on the brane. According to (2.5) we can accomplish this 
requirement by setting the boundary condition for the Bessel function: 

^2(r)|,=o = 0- (2-7) 

Since J2 is an oscillatory function, in the case of increasing (decreasing) warp factor the 
function (2.6) can have finite (infinite) number of zeros. Thus the above boundary condition 
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can be written in the form which quantizes the oscillation frequency w of the standing wave 
in terms of the curvature scale a, i.e. 



V-,=Xn, (2.8) 
\a\ 

where X„ is the n*'* zero of J2{f)- Correspondingly, the nodes of the standing wave in the 
bulk, the points where the functions (p{t, r) and u{t, r) vanish, can be considered as 4D 
space-time 'islands', where the matter particles are assumed to be bound. 

We consider the case with the increasing warp factor (a > 0) when the function 
J2(e~'^l'"l) in (2.5) has a finite number of zeros. For simplicity, in this paper we explore the 
case when J2 has a single zero at the position of the brane, r = 0, i.e. we assume that 

-=Xi^ 5.14 . (2.9) 
a 

In the equations of matter fields the oscillatory function (2.5) enters via some expo- 
nential functions: 



e 

^ — ' ni 

n=0 

where 6 is a constant. We suppose that the frequency ui of standing waves in the oscillatory 
metric function u{t, r) is much larger than the frequencies associated with the energies of 
particles on the brane. In this case we can perform the time averaging of oscillating 
exponents in the equation of the matter fields. 
Using the expression: 



00 



(m = 2n -h 1) 

sin™(^t)a!t = ^, ^ _ (2.11) 

27r 7 1^ 2'"(m/2)!^ - ^^^) 







we find that the time averages of some oscillatory functions are zero [28]: 

(^) = (u'> = (n) = ((e^")')=0, (2.12) 



where the primes and dots mean the derivatives with respect to the extra coordinate r and 
time coordinate t, respectively. We also obtain the simple formulas for the non- vanishing 
averages: 



^ 22"n!2 

n.=0 

(6Z)2"+i 



(6Z)2" 

n.=0 



6^' V -^T^ — = b^ZZ' \Iq + - (IqLi - IiLq)] , (2.13) 



22^^(271 + l)n!^ 

n=0 ^ ' 



where /q, Ii and Lq, Li are the modified Bessel and Struve functions of the argument bZ, 
respectively. 
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3 Localization of bosons 



Let us consider the localization problem of quantum fields with spin 0, 1 and 2 [28, 29]. 
The massless scalar field can be defined by the 5D action: 



5$ = y d'x^ g^'^dM^dN^ . (3.1) 



The corresponding Klein-Gordon equation, 

^ dM (v^5^^9iv$) = , (3.2) 

on the background (2.2) has the form: 

[d^ - e-"(52 + d^) - e^^'dl] $ = e^-^l"^! (e^^l^l^')' . (3.3) 
We look for the solution of this equation in the form: 

$(x^r) =e*?''^^'C(r) (3.4) 
(Greek letters are used for 4D indices), which transforms (3.3) into the equation, 

(e^'^Hc')' = e2»H M + vDe'^ + py^" - pj\ i . (3.5) 

When the frequency oj of standing waves is much larger than the frequencies associated 
with the energy pt of particles on the brane we can replace the oscillatory exponents in 
(3.5) by their time averages (2.13). Then we obtain the equation for the extra factor ^(r): 

which contains the r- dependent momentum: 

P\r) = ((e-) - 1) {pI+pD + ((e^") - l) p^ (3.7) 

On the brane, where « ~ 0, the parameters p^ can be regarded as components of 
energy-momentum along the brane which obeys the dispersion relation: 

pI-pI-pI-pI = ^ ■ (3.8) 

To study a general behavior of the extra part of the scalar zero mode wave function 
we explore (3.6) in two limiting regions: far from and close to the brane, respectively. The 
r -dependent momentum (3.7), which describes oscillatory properties of standing waves, has 
the following asymptotical forms: 

p2(r)| ~ const . (3.9) 
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Then we solve (3.6) in this hmiting regions and find: 



C{r)\r^o ~ const , 

Hence ^(r) has maximum on the brane and fahs off at the infinity as e"'^"'^'. 

In the action of scalar fields (3.1) the determinant (2.3) and the metric tensor with 
upper indices give the total exponential factor e^"'^', which obviously increases for a > 0. 
This is the reason why in the original brane models [3, 4] the scalar field zero modes with 
the constant r-dependent extra part can be localized on the brane only in the case of 
decreasing warp factor (i.e. a < 0). In our model the extra part of wave function (3.10) 
is not constant, moreover, for a > it contains the exponentially decreasing factor e"'*"'''!. 
For such an extra dimension factor the integral over r in the action (3.1) is convergent, 
hence 4D scalar fields are localized on the brane [28]. 

It is known that the transverse traceless graviton modes obey the equation of a massless 
scalar field in a curved background. Indeed, let us consider the metric fiuctuations: 

ds'^ = e^'^l'^l {g^^ + V) dx^'dx" - dr"^ , (3.11) 

where g^u is the metric tensor of the 4D part of (2.2): 

5,.. = (l,-e^-e^-e-2") . (3.12) 

Close to the brane (u ~ 0) for time averages of oscillatory exponents (2.13) we can use 
the approximation: 

« l + (t,) . (3.13) 

Thus the functions (n) can be regarded as r-dependent additive terms of Then the 
equations of motion for the fiuctuations h^i,: 

Om (V^/'^^^ivV) = , (3.14) 

are equivalent to the equation of motion of a scalar field (3.2) if we replace $ with h^,^. 
Accordingly, the condition of localization of spin-2 graviton field might be equivalent to 
that of spin-0 scalar field considered above. 

A little bit more complicated task is to show the localization of vector field zero modes 
[29]. Let us consider only the U{1) vector field (the generalization to the case of non-Abelian 
gauge fields is straightforward). 

The 5D action of vector fields, 

Sa = -\J g^'''g''''FMpFNR , (3.15) 

where 

Fmp = OmAp - dpAM , (3.16) 
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leads to the system of five equations: 

^9M(V^/'V^FNii)=0. (3.17) 

We seek for the solution of the system (3.17) in the form: 

ylt(x^) = v{r) atixn , 
A,{x^') = e"(*''^)i;(r) a,{x'') , 

Ay{x^) = e"(*''^)t'(r) ayix") , (3.18) 
A,{x^) = e-2"(*'^)i;(r) a^(x^) , 
Arix^) = , 

where a^{x'^) denote the components of 4D vector potential and scalar factor v{r) depends 
only on the extra coordinate r. 

Taking into account the equalities (2.12), the time averaging of the fifth equation of 
the system (3.17), 

= , (3.19) 

yields the Lorenz-like gauge condition: 



5" 



^do^Ap = rff'do.afi = , (3.20) 



where 7^^^ denotes the metric of 4D Minkowski space-time. The equation (3.20), together 
with the last expression of (3.18), can be considered as the full set of imposed gauge 
conditions. 

The remaining four equations of the system (3.17), 

d, (<7^^/"F5«) -^{Va Z"^.)' = , (3.21) 
after the time averaging and the use of (2.12) and (3.20) reduce to: 

V g^'^dadgap + e-^'^l"! (e^'^l'lt;')' = . (3.22) 

We require the existence of 4D vector waves localized on the brane, 

{xn ~ e^e'P-^^'' , (3.23) 

where pi, are the components of energy-momentum along the brane. Then the system 
(3.22) yields the single equation for v{r): 

(^e'^'^l^^v'y - P'^{r)v = , (3.24) 

where P'^{r) is done by (3.7). The solutions of (3.24) close to and far from the brane are: 

^j{r)\r^o ~ const , 

v{r)l^oo ~ e-^'^H . (3.25) 
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We see that the extra factor v{r) of the vector field zero mode wave function has a maximum 
on the brane and fahs off at the infinity as e~^''l''L 

In 5D brane models the implementation of pure gravitational trapping mechanism of 
vector field particles remains the most problematic. The reason is that in the vector action 
(3.15) the extra dimension parts of the determinant (2.3) and two metric tensors with upper 
indices cancel each other. Because of this in the original brane models the vector field zero 
modes (with the constant extra dimension part) cannot be localized on the brane for any 
sign of a. In our model the extra dimension part of vector field v{r) is not constant; at the 
infinity it decreases exponentially (3.25). Thus the integral over r in the action (3.15) is 
convergent, what means that 4D vector fields are localized on the brane [29]. 

4 Localization of fermions 

Now we investigate the localization problem for massless fermions [30]. 

For Minkowskian 4x4 gamma matrices {{'y'^,J^} = 2r]°^^) we use the Weyl basis: 





7* = ; . , 7' = : , (4.1) 



where / and a* [i = x, y, z) denote the standard 2x2 unit and Pauli matrices respectively. 
5D gamma matrices T"^ = h^T^ can be chosen as: 

ry = e-'^l'^l-"/2 , (4.2) 

_ ^-a\r\+u Y 

V = i-f^ . 

The 5D Dirac action for massless fermions: 

S^ = J d'-x^ [x^) T^'Dm^ [x^) , (4.3) 
contains the covariant derivatives: 

DA = dA + . (4.4) 

The non-vanishing components of the spin-connection in the background (2.2) are: 

f^f = - (e^l'^l-")' , (4.5) 
af = = ( e«/2' 

nf = (e--) 
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The corresponding to (4.3) 5D Dirac equation reads: 

iT^D^^' = i (T^Z)^ + T'Dr) ^ = . (4.6) 

For wave function of the bulk fermion field we use the chiral decomposition: 

^ (x^ r) = iPl {x") A(r) + iPr (x^) p{r) , (4.7) 

where A(r) and p(r) are the extra dimension factors of the left and right fermion wave 
functions respectively. We assume that 4D left and right Dirac spinors: 

7^Vl = -V'L , 7Vi? = +4^R , (4.8) 

correspond to the zero mode wave functions, i.e. they satisfy the free Dirac equations: 

il'^df.ijL = il^d^i^R = . (4.9) 

The solutions of (4.9) in our representation (4.1) can be written in the form: 



R 



V'l(x^)= (4.10) 

where the constant 2-spinors L and R satisfy: 

{pt + a'pi)L={pt-a'pi)R = ^ . (4.11) 

When the frequency oj of standing waves is much larger than the frequencies associated 
with the energies pt of the fermions on the brane we can time average the oscillatory 
functions in the Dirac equation (4.6). Time averages of the Dirac operators are: 

(iT'Dt) = ie-^l^^l 7*9t - sgn{r)-t\ 

{iT^D,) = ie-'^l'^l (e") 7^5, - sgn{r)^^ (4.12) 

and the equation (4.6) takes the form: 

i Wdt + (e"/2\ (^^'a^ + ^Vdy) + (e-"> Ydz] = e'''^^^ [2a sgn{r) + 9^] ^ . (4.13) 



Using the solutions of free equations (4.10) and the relations (4.11) it can be rewritten as 
the system: 



-e^l^'l [2a sgn(r) + dr\ cf'Vi{r) \ ( p{r)R 

-a^Viir) e^l^l [2a sgn{r) + d^] ) \ X{r)L 



. (4.14) 
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Here we have introduced the functions Vi{r): 

r.ir) = ((e-"/2\ _ i) = [I, (|C|Z/2) - l]p,, 



l)Py=[l0{\C\Z/2)-l]py, 



V.{r) = {{e-)-l)p, = [Io{\C\Z)-l]p, 



(4.15) 



where Z{r) is defined in (2.6). These functions, as (3.7), can be considered as the compo- 
nents of 'r -dependent momentum' of the spinor field: 

r'^{r) = r^ + r^ + r^ . (4.16) 

From the second equation of the system (4.14) it is straightforward to find 

P(^)^ = e"''' [2« sgnir) + dr] X{r)L . (4.17) 

Inserting (4.17) into the first equation of (4.14) and multiplying the result by a^Vi, we 
receive the second order differential equation for the function A(r): 



5a sgn{r) - — 



A' + 



4a(5(r) + 60^ - 2a sgn{r)— - p2g"2a|r| 



A = 



(4.18) 



Now, as for the case of bosonic fields, let us investigate this equation in the domains far 
from and close to the brane, respectively. 

Close to the brane the 'r -dependent momentum' (4.16) behaves as: 

V[r)l^^^ = Ar^ + 0{r^) , (4.19) 

where A is a constant, and the equation (4.18) takes the following asymptotic form: 



A" + 



5a sgn{r) 



A' + 



2 4a 

4ao(r) + 6a sgnyr) 



A = . 



The main solution of this equation in our approximation is: 

A(r)|,^±o = i?e-2»M , 



(4.20) 



(4.21) 



where i? is a constant. 

Note that, as a consequence of (4.17) and (4.21), in our setup the right fermionic modes 
are absent on the brane: 

/o(r-)|r-.±o = . (4.22) 



In the second limited region - far away from the brane, 

V{r)\r^±oo ~ const , 
and the equation (4.18) takes the asymptotic form: 

A" + 5a sgn{r)X' + 6a^X = , 



(4.23) 



(4.24) 
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with the solution: 

A(r)|._>±oo~e-3-IH . (4.25) 

Using (4.25) from the relation (4.17) we find also the asymptotic behavior of the extra 
dimension factor of the right fermion wave function: 

p(r)|,^±oo~e-2-M . (4.26) 

So in our model the extra dimension part of the left spinor wave function (4.7) has the 
maximum at the origin, 

X{r)\r=o = B, (4.27) 

decreases from the brane, and turns into the asymptotic form (4.25) at the infinity. When 
r — > oo the determinant (2.3) in the action integral (4.3) increases as e'*'^l''L However, 
fiinfbein in (4.12) contribute e""'^' and the extra dimension factor of left fermions (4.25) 
is proportional to e~^'^l''L Thus the overall r-dependent part decreases as e-^'^l'^l, i.e. the 
integral over r in (4.3) is convergent and zero modes of left fermions are localized on the 
brane. 

Due to (4.26) the extra dimension part of right fermions cancels the factor from the 
determinant (2.3) in the action (4.3). The r-component of the Dirac operator (4.12) does 
not contain an extra decreasing factor from fiinfbein and the integral over r of the term 
i^i T^dr^ in (4.3) diverges, i.e. the zero mode wavefunctions of right fermions actually are 
not normalizable. 



5 Classical particles 

Finally we consider the motion of a classical particle, or a photon, which obey 5D geodesic 
equation of motion: 

where k is the parameter of trajectory. Non-zero components of 5D Cristoffel symbols for 
the metric (2.2) are: 



1 / o^.\ - — + — ^ 1 

2 



= n\r\' r* =-(p~'^'^\ r* = r* =-(p^\ 

tr "'\' \ 1 zz o V ) 1 XX yy r) V-- ) 1 



tx 


- 

-^ty 


1 

= 2^' 




pX 

rx 




J- tt 


4( 




r =-1 

^ '■zz 2 







\' + \u', Tl, = a\r\' -u', (5.2) 

pr _ pr _ _t. (p2a|r|+«\ 
XX yy 2 V / ' 

For simplicity we shall consider the motion in (rj;)-plane, i.e. we suppose 

dy = dz = 0. (5.3) 
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Then y and z are not dynamical variables and the system (5.1) consists of three independent 
equations: 



dk"^ ''dkdk 



\dk 



, 



d'^x dt dx dr dx 

d^ ^ ""dkdk ^ dA? " ' 



(5.4) 



dk^ 







The second equation of this system, after dividing by dx/dk, can be rewritten as: 



d_ 

dk 



It's first integral is: 



, I dx\ u dt / I ,/ u'\ dr 
^"'d^J + 2d^+l""^" +tJ dk^ 



^ - T/p-«k-|-"/2 



, 



(5.5) 



(5.6) 



where the constant V corresponds to the component of the particles velocity along the 
brane. 

Inserting (5.6) into the first equation of the system (5.4) we find: 



dk 



2 dt 



(5.7) 



In the case of fast oscillations of standing waves the time average of the last term yields 
zero (since {u) = 0) and 

-'^I^l . (5.8) 



dt 
dk 



The integration constant in this expression was included in the definition of k, so that on 
the brane (r = 0) the parameter k coincides with the coordinate time t. 

Inserting (5.2), (5.6) and (5.8) into the last equation of (5.4) and multiplying it by 
dr/dk we find: 

d 



dk 







(5.9) 



Time average of the last term in this equation gives zero {{u) = 0) and the first integral is: 



1 /dr 

2 (dk 



+ a{l-V^)\r\ = e , 



(5.10) 



where the constant e > corresponds to the energy per unit mass and a(l — y^)|r| plays 
the role of the trapping gravitational potential. 

To find a connection of the parameter k with the proper time let us insert (5.3), (5.6), 
(5.8) and (5.10) into the definition of the interval (2.2), 

ds^ = e^^l^'ldi^ - e^'^l'^l+^dx^ - dr^ = [(1 + 2a|r|)(l - V^) - e] dk^ . (5.11) 

We see that on the brane (r = 0) for photons e = 1 and for massive particles < e < 1. 
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From (5.10) it is clear that the motion towards the extra dimension r is possible when 

e-a(l-y2)|r| >0 , (5.12) 

and for any energy e there exists the maximal distance in the bulk, 

\r\max ~ - , (5.13) 
a 

the particle can reach, i.e. the classical particles are trapped on the brane. 

In the standard brane approach with decreasing warp factor (a < 0) localization was 
achieved due to the fact that the extra space actually is finite [3, 4]. In our case the 
increasing of the brane warp factor, e^"!'^ ', creates the potential well that confines particles. 

Another point is that we had neglected the influence of the oscillating exponents in 
(2.2). In this approximation the anisotropy of the background metric for classical particles 
is hidden. However, anisotropic nature of the primordial metric (2.2) can be exhibited in 
cosmological solutions [31, 32]. 



6 The summary and discussions 

In this letter we have demonstrated the pure gravitational localization of all kinds of matter 
fields within the 5D standing wave braneworld [13, 14]. The main differences of our model 
from the standard brane approaches [3, 4] are: 

i) Metric ansatz contains the increasing warp factor; 

ii) In the bulk there exist rapidly oscillating standing waves which localize fields on 
the brane; 

iii) The extra dimension factors of zero modes are not constant and far from the brane 
behave as: 



^>(r) 


~ e~ 


-4a|r 




~ e~ 


-ia\r 




~ e~ 


-3a r 




~ e~ 


-2a\r 


hMN{r) 


~ e~ 


-ia\r 



where A , ^l, ^'/j and Kmn correspond to the scalar, vector, left and right fermion 
and graviton respectively. 

The advantage of our model is that the localization mechanism is universal, i.e. it 
should work also for interacting fields. If we suppose an extra dimensional profile for 
interacting fields as displayed in (6.1), then it is clear that the integrals over r of standard 
5D interacting terms like ^A^jj\f, y^<I>^/j^2., ^/g^'^ A]\f , etc., are convergent and the 
extra dimensional factors only lead to renormalization of 4D coupling constants. 
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